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I discuss mechanisms of heavy quark production in (real) photon-nucleon 
and (real) photon - (real) photon collisions. In particular, I focuse on ap- 

■ plication of the Saturation Model. In addition to the main dipole-nucleon 
\Q ' or dipole-dipole contribution included in recent analyses, I propose how 

, to calculate within the same formalism the hadronic single-resolved con- 

tribution to heavy quark production. At high photon-photon energies this 
yields a sizeable correction of about 30-40 % for inclusive charm production 

■ and 15-20 % for bottom production. Adding all possible contributions to 
O |' e + e~ —* bbX together removes a huge deficit observed in earlier works but 
^ ■ does not solve the problem totally. 
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1. Introduction 

The total cross section for virtual photon - proton scattering in the re- 
gion of small x and intermediate Q 2 can be well described by the Saturation 
Model (SAT-MOD) 1 . The very good agreement with experimental data 
can be extended even to the region of rather small Q 2 by adjusting an effec- 
tive quark mass. At present there is no deep understanding of the fit value 
of the parameter as we do not understand in detail the confinement and the 
underlying nonperturbative effects related to large size QCD contributions. 

In this presentation I shall limit to the production of heavy quarks which 
is simpler and more transparent for real photons. Here one can partially 
avoid the problem of the poor understanding of the effective light quark 
mass, i.e. the domain of the large (transverse) size of the hadronic system 
emerging from the photon. 

It was shown recently that the simple SAT-MOD description can be 
succesfully extended also to the photon-photon scattering :2;. The heavy 
quark production in photon-photon collisions is interesting in the context 
of a deficit of standard QCD predictions relative to the experimental data 
as observed recently for b quark production. 
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2. Heavy quark production in photon-nucleon scattering 

In the so-called dipole picture the cross section for heavy quark-antiquark 
(QQ) photoproduction on the nucleon can be written as 

a lN ^ QQ (W) = J d 2 pdz \^ Q (p,z)\ 2 a dN (P,z,W) , (1) 

where <E>t is (transverse) quark-antiquark photon wave function (see for 
instance @J) and a^N ls the dipole- nucleon total cross section. Inspired by 
its phenomenological success ^ we shall use the SAT-MOD parametrization 
for (JdN- Because for real photoproduction the Bjorken-x is not defined we 
are forced to replace xhy x g 

In Fig. la I show predictions of SAT-MOD for charm photoproduction. 
The dotted line represents calculations based on Eq.(JJ). The result of this 
calculation exceeds considerably the fixed target experimental data. One 
should remember, however, that the simple formula Q applies at high en- 
ergies only. At lower energies one should include effects due to kinematical 
threshold. In the momentum representation this can be done by requiring: 
Mqq < W, where Mqq is the invariant mass of the final QQ system. This 
upper limit still exceeds the low energy experimental data. There are phase 
space limitations in the region x g — > 1 which have been neglected so far. 
Those can be estimated using naive counting rules. Such a procedure leads 
to a reasonable agreement with the fixed target experimental data. 

The deviation of the solid line from the dotted line gives an idea of 
the range of the safe applicability of SAT-MOD for the production of the 
charm quarks/antiquarks. The cross section for W > 20 GeV is practically 
independent of the approximate treatment of the threshold effects. SAT- 
MOD seems to slightly underestimate the HI collaboration data 5 . For 
comparison in Fig.l we show the result of similar calculations in the collinear 
approach (thick dash-dotted line) with details described in [3]. 

The calculation above is not complete. For real photons a vector domi- 
nance contribution due to photon fluctuation into vector mesons should be 
included on the top of the dipole contribution. In the present calculation 
we include only the dominant gluon-gluon fusion component. Then 

tf™Qn(W) = ^2 % [ dx v dx N g v (x v ,fj, F ) 9n(xn,Hf) (7 gg ->QQ(W) ■ 
v J v J 

(2) 

Here the fy constants describe the transition of the photon into vector 
mesons (p, u>, </>). The gluon distributions in vector mesons are taken as 
that for the pion [7j. 

The dash-dotted line in Fig. la shows the VDM contribution calculated 
in the leading order (LO) approximation for a gg _^QQ. The so-calculated 
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Fig. 1. The cross section for j + p -> QQ^- The dotted line: standard SAT-MOD, 
the dashed line: includes kinematical threshold, the solid line includes in addition 
a suppression by (1 — cc c ) 7 , the thick dash-dotted line: collinear approximation and 
the thin dash-dotted line: the LO VDM contribution. 



VDM contribution is small at small energies but cannot be neglected at 
high energies. 

The situation for bottom photoproduction seems similar. In Fig. lb I 
compare the SAT-MOD predictions with the data from the HI collaboration 
[Hj. Here the threshold effects survive up to very high energy W ~ 50 GeV. 
Again the predictions of SAT-MOD are slightly below the HI experimental 
data point. The relative magnitude of the VDM component is similar as for 
the charm production. 

3. Heavy quark production in photon-photon scattering 

In the dipole-dipole approach the total cross section for 77 — > QQ pro- 
duction can be expressed as 

o%^qq(W) = Y, I ^ Q0 (pi^i)W 2f2 (p2,Z2)\ 2 a dd (pi,P2,x Qf ) d 2 Pl dz 1 d 2 p 2 dz 2 
h+Q J 

+ Yl / \® flIl (Pii z i)\ 2 \® QQ (P2i z 2)\ 2(7 dd(pi,P2,XfQ) d 2 pxdzxd 2 p 2 dz 2 , 
h+Q J 

(3) 

where add is the dipole-dipole cross section. 



4 



EPIPH03 PRINTED ON FEBRUARY 1, 2008 



There are two problems associated with direct use of (jHJ). First of all, it 
is not completely clear how to generalize add from a d N parametrized in pQ. 
Secondly, formula (j2J) is correct only at W S> 2mg. At lower energies one 
should worry about proximity of the kinematical threshold. 

In a very recent paper [2] a new phenomenological parametrization for 
add has been proposed. The phenomenological threshold factor in [2j does 
not guarantee automatic vanishing of the cross section exactly below the 
true kinematical threshold W = 2m a + 2m;,. Therefore, instead of the 
phenomenological factor I impose an extra kinematical constraint: M^j + 
Mqq < W on the integration in ©. 

It is also not completely clear how to generalize the energy dependence of 
a dN in photon-nucleon scattering to the energy dependence of add in photon- 
photon scattering. In [Hj I have defined the parameter which controls the 
SAT-MOD energy dependence of add m a symmetric way with respect to 
both photons. In comparison to the prescription in |2], our prescription 
leads to a small reduction of the cross section far from the threshold [3] . 

Up to now we have calculated the contribution when photons fluctuate 
into quark-antiquark pairs, which is not complete. The dipole approach 
must be supplemented by the resolved photon contribution. This contribu- 
tion can be estimated in the vector dominance model when either of the 
photons fluctuates into vector mesons. If the first photon fluctuates into 
the vector mesons, the so-defined single-resolved contribution to the heavy 
quark-antiquark production can be calculated analogously to the photon- 
nucleon case as 

< 7 %q(W0 =EtT / \^ Q (P2,z 2 )\ 2 a Vld (p 2 ,x 1 ) <Pp 2 dz 2 , (4) 

where ay ld is vector meson - dipole total cross section. In the spirit of SAT- 
MOD, we parametrize the latter exactly as for the photon-nucleon case [Q 
with a simple rescaling of the normalization factor a$ v = 2/3 • a$ N . In 
the present calculation a^ N as well as the other parameters of SAT-MOD 
are taken from 1 . Analogously, if the second photon fluctuates into vector 
mesons we obtain 

^qq(W) = J2^- I \*F(pi,*)\ 3 **h(pi>**) (5) 

y 2 J V2 J 

This clearly doubles the first contribution @ to the total cross section. 

The integrations in @ and (JSJ) are not free of kinematical constraints. 
When calculating both single-resolved contributions, it should be checked 
additionally if the heavy quark-antiquark invariant mass Mqq is smaller 
than the total photon-photon energy W (see 0). 
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Fig. 2. Different contributions to the inclusive charm (left panel) and bottom (right 
panel) production. The long-dashed line: the dipole-dipole contribution, the dash- 
dotted line: the single-resolved contribution, the lower dashed line: the 2Q2Q 
contribution, the dotted line: the direct contribution, the gray solid line: double- 
resolved contribution. The experimental data for inclusive c/c production are from 
Ref.0. 



In Fig. 2 I show different contributions to the inclusive c/c (left panel) 
and b/b (right panel) production in photon-photon scattering. The thick 
solid line represents the sum of all contributions. 

Let us start from the discussion of the inclusive charm production. The 

experimental data of the L3 collaboration j^j are shown for comparison. 

The modifications discussed above lead to a small damping of the cross 

section in comparison to Ref.|5]. The corresponding result (long-dashed 

line) stays below the recent experimental data of the L3 collaboration 

The hadronic single-resolved contribution constitutes about 30 - 40 % of 

the main SM contribution. At high energies the cross section for the 2c2c 

component is about 8 % of that for the single cc pair component. In the 

inclusive cross section its contribution should be doubled because each of 

the heavy quarks/ antiquarks can be potentially identified experimentally. 
At higher energies the direct contribution is practically negligible. In 

contrast, the hadronic double-resolved contribution, when each of the two 
photons fluctuates into a vector meson jj! is shown by the thin solid line in 
the figure becomes important only at very high energies relevant for TESLA. 
Here we have consistently taken gy(xv,H 2 p) = g 1T {xv : ^ 2 F )- 

The situation for bottom production (see right panel) is somewhat dif- 
ferent. Here the main SAT-MOD component is dominant. Due to smaller 
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charge of the bottom quark than that for the charm quark, the direct com- 
ponent is effectively reduced with respect to the dominant SAT-MOD com- 
ponent by the corresponding ratio of quark/antiquark charges: (1/9) 2 : 
(4/9) 2 = 1/16. The same is true for the 2626 component. Here, in addition, 
there are threshold effects which play a role up to relatively high energy. 

3.1. Short- versus long-distance phenomena 

What are typical distances probed in heavy quark production ? Is the 
heavy quark production a short distance phenomenon ? These questions 
can be easily answered in the mixed representation formulation considered 
in the present paper. In Fig|3]I display the integrand of 

a 11 ^ C c(W) = J I{p\, pi) dpidpi . (6) 

The maxima in Fig |3] correspond to the most probable situations. For one 
light (m u = nid = mo, m s = ttiq + 0.15 GeV) and second heavy quark- 
antiquark pair the map is clearly asymmetric. One can observe a ridge 
parallel to the p\ or p2 axis. There is no well localized maximum. Both 
short and long distances are probed. 

For comparison, in the bottom part of the figure, I show similar maps 
when both pairs consist of light (u,d,s) quarks/antiquarks (left-bottom) 
and in the case when both pairs consist of charm quarks/antiquarks (right- 
bottom). For light quarks (u,d,s) one observes a clear maximum at (pi,p2) 
= (1 GeV -1 , 1 GeV -1 ) = (0.2 fm, 0.2 fm). In this case a non-negligible 
strength extends, however, up to large distances p\ and p2- Only in the case 
of the production of two cc pairs, the cross section is dominated exclusively 
by short-distance phenomena. 

3.2. e+e" -> e+e - 66V 

Up to now no attempt was done to unfold experimentally the cross 
section for 77 — > bbX. Only the cross section for the e + e~ — ► bbX reaction 
was obtained recently by the L3 and OPAL collaborations at LEP2 [51 ITU]. 
The measured, positron and electron antitagged cross sections cannot be 
described as a sum of direct and single-resolved contributions, even if next- 
to- leading order corrections are included dl]]. The measured cross section 
exceeds the theoretical predictions by a large factor. This is a new situation 
in comparison to the charm production where the deficit is much smaller. 

The cross section for the e + e~ — ► 66 X reaction when both positron and 
electron are antitagged can be easily estimated in the equivalent photon 
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Fig. 3. A map of ^^Ypldfz'^ at W = 100 GcV for thc first ( lcft - to P P ancl ) and 
the second (right-top panel) photon fluctuating into cc. For comparison I show 
analogous maps for light quark-antiquark pairs (left-bottom panel) and for the 
case when both pairs consist of charm quarks/ antiquarks (right-bottom panel). 



approximation (EPA) as 

a(e + e~ ->■ bbX;W ee ) = J dx A dx B fA(E b ,6 max ,x A )fB(Eb,0 max j %B ) 

a(77 bbX; W) , (7) 

where f A and /# are virtuality-integrated flux factors of photons in the 
positron and electron, respectively, and 9 max is the maximal angle of the 
positron/electron not, to be identified by the experimental apparatus. In 
the present analysis we calculate the integrated flux factors f A and fs in 
a simple logarithmic approximation. The photon-photon energy can be 
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Fig. 4. The dependence of the integrand of Eq.Q on the photon-photon energy 
W 77 for xf — (solid) and xp = ± 0.5 (dashed) for the bb production for the 
direct mechanism (left panel) and in the dipole-dipole scattering in the Saturation 
Model (right panel) with the present prescription for the energy dependence of the 
dipole-dipole cross section. In this calculation W ee — 190 GeV. 



calculated in terms of photon longitudinal momentum fractions xa and 
xb in the positron and electron, respectively, as W = ^XAXBS ee - It is 
instructive to visualize how different regions of W 77 contribute to a(e + e~ — > 
bbX;W ee ). For this purpose it is useful to transform variables from xa,xb 
to xp = xa — xb and W T7 = W. Then 

a(e + e~ — ► bbX; W ee ) = J dxpdW J /a(-^6, G max , XA)fs(Eb, "max j Xb ) 

ct(77 bbX;W) , (8) 

where the Jacobian J is a simple function of kinematical variables. 

The integrand I(xf,W) of (jHJ) for bb is shown in Fig|l] for the direct 
production (left panel) and for the Saturation Model (right panel) including 
all contributions considered in the present analysis. Quite a different pattern 
can be observed for the two mechanisms. While for the direct production one 
is sensitive mainly to low-energy photon-photon collisions, in the Saturation 
Model the contributions of high-energies cannot be neglected and one has 
to integrate over W 77 essentially up to W ee . This difference in W 77 is 
due to different energy dependence of 17(77 ~~ > bb; W^y) for the different 
mechanisms considered, as has been discussed above. Even in the latter 
case the integrated cross section is very sensitive to the region of not too 
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Table 1. Cross sections in pb for e + e — ► bbX for LEP2 averaged energy W ee = 
190 GeV. 



direct 


bb 

SAT-MOD 


2b2b 
SAT-MOD 


SR 
SAT-MOD 


sum 


L3 


OPAL 


1.21 


6.1-7.4 


0.034 


1.92 


9.3-10.6 


13.1 ± 2.0 ± 2.4 


14.2 ± 2.5 ± 5.0 



high energies W ~ 20 GeV, where the not-fully-understood threshold effects 
may play essential role. 

For LEP2 averaged energy < W ee >~ 190 GeV the cross section inte- 
grated taking into account experimental cuts is a(e + e~ — ► bbX) = 6.1 pb 
(C = 1) or 7.4 pb (C = 1/2) 3 for the dipole-dipole SAT-MOD scatter- 
ing process. The corresponding cross section for the direct production is 
a(e + e~ — * bbX) = 1.2 pb. The hadronic single-resolved contribution calcu- 
lated here in the Saturation Model is very similar in size to that calculated 
in the standard collinear approach ^Uj. As can be seen in Table 1 the 
2b2b contribution is practically negligible. We have completely omitted the 
double-resolved contribution which is practically negligible. The sum of the 
direct, bb SAT-MOD, 2626 SAT-MOD and the hadronic single resolved SAT- 
MOD component is 9.3-10.6 pb in the case when no transverse momenta 
cuts on the main SAT-MOD component are included and 6.4-7.1 pb with 
the cuts. These numbers should be compared to experimentally measured 
a(e + e~ -» bbX) = 13.1 ± 2.0 (stat) ± 2.4 (syst) pb QI] (L3) and prelimi- 
nary a{e + e~ -> bbX) = 14.2 ± 2.5 (stat) ± 5.0 (syst) pb [H] (OPAL). In 
comparison to earlier calculations in the literature, the theoretical deficit is 
much smaller. The success of the present calculation relies on the inclusion 
of a few mechanisms neglected so far - in particular the dipole-dipole con- 
tribution which, in our opinion, is not contained in the standard collinear 
approach. 

Up to now only the W 77 -integrated cross section has been determined 
experimentally. This, in fact, does not allow to identify experimentally 
whether the problem is in low or high W 77 . In order to identify better 
the region where the standard collinear approach fails it would be useful to 
bin the experimental cross section in the intervals of W 77 making use of a 
possibility to measure W V i S which can be related to W 77 via a suitable Monte 
Carlo program. At present, even splitting the cross section for e + e~ — ► 66A 
into cr(W 77 < Wo) and a(W 11 > W ) for W ~ 20 GeV would be useful 
and should shed more light on the problem of the experimental excess of 66 
relative to the "standard" QCD approach. 
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3.3. Quark- antiquark correlations 

So far mainly the integrated cross section for heavy quark/antiquark 
production was considered in the literature. Only in a few cases inclusive 
distributions in transverse momentum or rapidity (see e.g. |12j) were pre- 
sented. No attempts have so far been made to analyze the final state in 
more detail. In our opinion investigating correlations between heavy quark 
- heavy antiquark could be much more conclusive in identifying the produc- 
tion mechanisms than the integrated cross section or even a single variable 
distribution. 

In principle, any correlation between two kinematical variables of the 
final quark and antiquark would be of interest. We suggest that the following 
final quark/antiquark momentum fractions: 

_ PQ > 
Q ~ \p Q \ 71 ' 

X Q ~ i^TT n "n ' v y J 

\pq\ 

where 

*» = J*T (10) 

would be very useful to separate the different mechanisms (approaches). 
In the definition above pq and pq are momenta of the heavy quark and 
antiquark, respectively, and p 71 is the momentum of the first photon, all 
in the photon-photon center-of-mass frame. By definition -1 < xq,xq < 
1. Similar quantities are being used at present when analyzing, e.g., jet 
production at HERA to separate out resolved and direct processes. 



In Fig|5] I present a sketch of naive expectations. Although a precise 
map requires detailed calculations for each mechanism separately, which is 
beyond the scope of the present analysis, it is obvious that the separation 
of different mechanisms here should be much better than for any inclusive 
spectra. In the case of dipole-dipole approach (the elongated ellipses) this 
would require going to the momentum representation. The mixed represen- 
tation used in the present paper is useful only for integrated cross sections. 

Experimentally, the analysis suggested would be difficult at LEP2 be- 
cause of rather limited statistics. We hope that such an analysis will be 
possible at the photon-photon option at TESLA. At present, even localiz- 
ing a few LEP2 coincidence cc events in the diagram x c versus x s would be 
instructive. 




Fig. 5. The expected locii in (xq,xq) space of different mechanisms considered in 
the present analysis. SRi / SR2 means that the first/second photon was trans- 
formed into vector mesons and DR means that each of the photons was transformed 
into a vector meson. The dashed circle is the locus corresponding to the pairs emit- 
ted from the middle of the gluonic ladder (not discussed in the text). 



4. Conclusions 

There is no common consensus in the literature on detailed understand- 
ing of the dynamics of photon- nucleon and photon-photon collisions. In this 
presentation I have limited the discussion to the production of heavy quarks 
simultaneously in photon-nucleon and photon-photon collisions at high en- 
ergies. The sizeable mass of charm or bottom quarks sets a natural low 
energy limit on a naive application of SAT-MOD. Here a careful treatment 
of the kinematical threshold is required. 

We have started the analysis from (real) photon-nucleon scattering, 
which is very close to the domain of SAT-MOD as formulated in pQ. If 
the kinematical threshold corrections are included, SAT-MOD gives simi- 
lar results as the standard collinear approach for both charm and bottom 
production. We have estimated the VDM contribution to the heavy quark 
production. 
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The second part of the present analysis has been devoted to real photon 
- real photon collisions. For the first time in the literature we have estimated 
the cross section for the production of 2c2c final state. We have found that 
this component constitutes up to 10-15 % of the inclusive charm production 
at high energies and is negligible for the bottom production. We have shown 
how to generalize SM to the case when one of the photons fluctuates into 
light vector mesons. It was found that this component yields a significant 
correction of about 30-40 % for inclusive charm production and 15-20 % for 
bottom production. We have shown that the double resolved component, 
when both photons fluctuate into light vector mesons, is nonnegligible only 
at very high energies, both for the charm and bottom production. 

I have shown that the production of cc pairs (the same is true for bb) is 
not completely of perturbative character and involves both short- and large- 
size contributions. The latter as nonperturbative are unavoidably subjected 
to some modeling. Present experimental statistics do not allow extraction of 
cross sections for the 77 — ► bb reaction and therefore it is not clear where the 
observed deficit resides. It is not excluded that the apparent deficit of bot- 
tom quarks may reside at photon-photon energies close to threshold. This 
is a region where the underlying physics has never carefully been studied. 

Finally I have discussed a possibility to distinguish experimentally the 
different mechanisms discussed in the present paper by measuring heavy 
quark - antiquark correlations. This suggestion requires, however, further 
detailed studies of the Monte Carlo type, including experimental possibilities 
and limitations. 

The present calculation is not fully consistent as far as parameters are 
considered. I have taken SAT-MOD parameters from Ref.pQ, where only 
dipole (quark-antiquark) component was fitted to the total cross section for 
7*p scattering and added the resolved photon in addition. In a consistent 
treatment one should include in addition the resolved photon component 
in the fit to the HERA structure function data. Such an analysis has been 
completed only recently [PS] . 
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